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Abstract 

Early life history stages of marine organisms are generally thought to be more sensitive to environmental stress than adults. 
Although most marine invertebrates are broadcast spawners, some species are brooders and/or protect their embryos in 
egg or capsules. Brooding and encapsulation strategies are typically assumed to confer greater safety and protection to 
embryos, although little is known about the physico-chemical conditions within egg capsules. In the context of ocean 
acidification, the protective role of encapsulation remains to be investigated. To address this issue, we conducted 
experiments on the gastropod Crepidula fornicata. This species broods its embryos within capsules located under the female 
and veliger larvae are released directly into the water column. C. fornicata adults were reared at the current level of C0 2 
partial pressure (pC0 2 ) (390 uatm) and at elevated levels (750 and 1400 uatm) before and after fertilization and until larval 
release, such that larval development occurred entirely at a given pC0 2 . The pC0 2 effects on shell morphology, the 
frequency of abnormalities and mineralization level were investigated on released larvae. Shell length decreased by 6% and 
shell surface area by 11% at elevated pC0 2 (1400 uatm). The percentage of abnormalities was 1.5- to 4-fold higher at 
750 uatm and 1400 |iatm pC0 2 , respectively, than at 390 |iatm. The intensity of birefringence, used as a proxy for the 
mineralization level of the larval shell, also decreased with increasing pC0 2 . These negative results are likely explained by 
increased intracapsular acidosis due to elevated pC0 2 in extracapsular seawater. The encapsulation of C. fornicata embryos 
did not protect them against the deleterious effects of a predicted pC0 2 increase. Nevertheless, C. fornicata larvae seemed 
less affected than other mollusk species. Further studies are needed to identify the critical points of the life cycle in this 
species in light of future ocean acidification. 
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Introduction 

Early life history stages of marine species, including embryos 
and larvae, are of crucial importance in population dynamics as 
they ensure dispersion, colonize new areas and sustain populations 
[1]. Their success in development and final recruitment are 
essential for the persistence of viable populations. Early stages of 
marine invertebrates are in general morphologically and ecolog- 
ically distinct from the adult stage and are generally thought to be 
more sensitive to environmental stress [1] although, in some cases, 
they may be more tolerant than adults, e.g. some Antarctic species 
exposed to warming [2]. In the context of climate change, early 
development may be affected by various factors, such as 
temperature increases, hypoxia zones or ocean acidification. Due 
to the increase in atmospheric pC0 2 predicted for the end of the 
century (from 475 to 1313 (Xatm according to the Intergoverne- 
mental Panel on Climate Change (IPCC)), pH in surface seawaters 
is likely to decline by 0.06 - 0.32 units [3], leading to a decrease in 
carbonate ion concentrations (C0 3 z_ ) and a reduction in the 
calcium carbonate saturation state (£1) [4] . Due to these changes in 
seawater carbonate chemistry, ocean acidification is considered a 



major threat to calcifying marine species, affecting their physiology 
and impairing their ability to build calcium carbonate shells and 
skeletons [8,9,10,11], which can ultimately modify their behavior 
and distribution [5,6,7]. Early life stages (embryos, larvae and 
juveniles) of calcifying species are thus expected to be highly 
affected by ocean acidification [12,13], as opposed to non-calcified 
larvae which are predicted to be more tolerant [14,15]. This 
relatively higher vulnerability is likely due to fragile larval skeletons 
[9] and their high ratio of exposed surface-to-body mass compared 
to adults [16]. Identifying life history stages that are the most 
vulnerable to global change is needed to determine bottlenecks for 
species persistence and addressing their sensitivity to acidification 
is a major issue in a changing ocean [9]. 

Responses to near-future (end of century) levels of/>C0 2 depend 
on species, populations, habitats and developmental stages 
[9,17,18,19] and understanding these effects on the early life 
stages requires taking into account the complete developmental 
cycle, from egg to juvenile [12]. In particular, the impact of 
elevated />CO z /decreased pH on early life stages has been 
investigated in a broad range of species, including corals [20,21], 
echinoderms [22,23,24], crustaceans [15,16,25], mollusks [26,27], 
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and fish [28,29]. In mollusks, which have been studied intensively 
(see refs. [26,30] for a review), deleterious effects of increased 
pG0 2 have been demonstrated on fertilization success [31,32], 
hatching success [33,34], larval survival [35,36], growth [37,38], 
shell formation [39,40], development duration [41,42], and 
settlement [43,44]. 

Most of the species studied are broadcast spawners, which may 
be considered particularly vulnerable to ocean acidification 
because fertilization and complete pelagic larval life occur in the 
water column [9,45]. Whether alternative reproductive modes are 
affected in a changing ocean is still poorly documented. Brooding 
and/ or egg laying in egg masses or capsules are typically assumed 
to confer protection to developing embryos [46,47]. For example, 
it has been shown that encapsulated embryos of some gastropod 
species survive better in conditions of salinity stress than embryos 
removed from their capsule [48,49]. A few studies have explored 
the effects of decreased pH on embryos brooded and/ or laid in 
benthic gelatinous egg masses or in egg capsules in bivalve [50], 
gastropod [35,41,51,52] or cephalopod [27,53] mollusks. Depend- 
ing on the study, reduced pH has different effects that are related 
to the range of species habitats and the strategy to protect embryos 
(brooding, egg masses, capsules) studied as well as the source of pH 
change (J>G0 2 increase, salinity stress). Encapsulation has been 
suggested to protect embryos against ocean acidification [41,54], 
whereby the buffer capacity of intracapsular fluids may reduce the 
potential effect of extracapsular elevated pC0 2 seawater. 

To study this issue in a non-broadcast-spawner species, we 
chose the slipper limpet Crepidula fomicata, Linne 1758 (Gastro- 
poda) as our biological model. Native to the northeast American 
coast, this species was introduced in Europe at the end of the 1 9 th 
century, primarily via oyster farming [55], and has now become 
invasive in bays and estuaries where it reaches very high densities 
of up to several thousands of individuals per m [56]. It has a 
bentho-pelagic life cycle, with a number of original features. 
Benthic adults form stacks with males at the top and females at the 
bottom. After internal fertilization, females brood their embryos in 
egg capsules for 3 to 4 weeks [57,58]. Capsules are protected 
between the neck and the propodium of the female parent and 
attached to the substratum to which the female is fixed [58]. Each 
female spawns between 28 and 64 capsules, each containing 300 to 
500 embryos [57]. At the end of capsular development, veliger 
larvae are released at a size of about 400 urn in length into the 
water column where they spend between 2 and 7 weeks [59,60]. 
Upon reaching competence (800-1000 urn in length), larvae are 
able to metamorphose and settle on hard substrata [59,60]. 

The objective of this work was to investigate the effects of near- 
future levels of pC0 2 on the development of C. fomicata 
encapsulated embryos by studying the shell morphology and 
mineralization level of released larvae. To ensure that the 
complete development, from the egg to the released larva, 
occurred under high pG0 2 , parents were conditioned to the 
different pC0 2 levels before mating occurred. 

Methods 

Crepidula fomicata adult collection and culture 

C. fomicata stacks were collected by SCUBA divers on 30 
November 2011, after the end of the reproductive period [61] in 
Morlaix Bay (northwestern Brittany, France), at the "Barre des 
Flots" site (3°53.015'W; 48°40.015'N). No specific permissions 
were required for sampling at the selected location, as it is not 
privately-owned or protected. Field sampling did not involve 
endangered or protected species. 



After being held 6 weeks in natural ambient unfiltered seawater, 
C. fomicata adults were randomly distributed into 18 aquarium 
tanks of 10 L each (adapted from [62]) and reared for 24 weeks (12 
January 2012 to 28 June 2012) in three pC0 2 treatments selected 
according to the recommendations of Barry et al. [63]: (1) 
390 uatm (pH on the total scale (pHx) = 8.07) as the current 
pC0 2 (control), (2) 750 uatm (pH T = 7.82) and (3) 1400 uatm 
(pH x = 7.56); the former two pC0 2 levels are pessimistic 
scenarios predicted for the end of the century by the IPCC [3]. 
The/>C0 2 was adjusted by bubbling C0 2 -free air (current control 
pC0 2 ) or pure C0 2 (two elevated pC0 2 treatments) in three 100 L 
header tanks supplied with unfiltered seawater pumped at the foot 
of the Station Biologique de Roscoff. Each of the three pC0 2 
treatments had six 10 L replicate aquaria. This was an open 
system, and C0 2 -treated seawater from the mixing header tanks 
was continuously supplied to the 18 aquaria (6 per pC0 2 
condition), at a rate of 9 L h -1 (i.e. a renewal rate of 90% h -1 ). 
Aquaria were placed in a thermostatic bath where temperature 
was controlled to within ±0.2°C using 150 to 250 W submersible 
heaters. C. fomicata adults were grown at four successive 
temperature levels (10, 13, 16 and 19°C) which corresponded to 
the range of in situ temperatures typically encountered in our study 
area (Service d'Observation de la Mer et du LITtoral data). Adults 
were reared for four weeks at each temperature level. Changes in 
temperature were implemented slowly, with increases of 0.2°C 
day 1 over a period of two weeks. 

pC0 2 and temperature were monitored and controlled by an 
off-line feedback system (IKS Aquastar, Karlsbad, Germany) that 
regulated the addition of gas in the header tanks and the on/ off 
heater switch in the thermostatic bath. The pH values of the 
system were adjusted from daily measurements of pHx in each of 
the 18 aquaria using a pH meter (826 pH mobile, Metrohm AG, 
Herisau, Switzerland) calibrated using Tris/HCl and 2-aminopyr- 
idine/HCl buffers [64]. Slipper limpets were fed three times a 
week with a mix made from a stock solution of Chaetoceros gracilis 
(~15xl0 6 cells mL ') and Isochrysis affinis galbana (~26xl0 6 cells 
mL~ ). This algal mix (400 mL) was distributed in each aquarium. 
Seawater flow was stopped for two hours to allow the limpets to 
feed. 

Seawater parameters were monitored throughout the experi- 
ment in each of the 18 aquaria. pHx and temperature were 
recorded daily. Total alkalinity was measured every four weeks, by 
0.01 N HC1 potentiometric titration on an automatic titrator 
(Titroline alpha, Schott SI Analytics, Mainz, Germany) in a 
20 mL seawater sample taken from each aquarium. Salinity was 
also measured every four weeks with a conductimeter (LF 330/ 
SET, WTW, Weilheim, Germany) and it varied between 
34.2 ±0. 1 and 35.1±0.1 over the course of the experiment. The 
carbonate chemistry of seawater, i.e. dissolved inorganic carbon 
(DIC), exact pC0 2 and saturation state of aragonite (Aat) were 
calculated for each pC0 2 and temperature treatment using 
C0 2 SYS software [65] with the constants of Mehrbach et al. 
[66] refitted by Dickson and Millero [67]. Mean values of these 
parameters are given in Table 1. 

Larvae collection 

In C. fomicata, at the end of the embryonic development, the 
capsule membrane splits and veliger larvae are released in 
seawater. To prevent released larvae from escaping from their 
source aquarium, 200 urn mesh size nets covered the overflow 
outlet of each aquarium. Offspring presence was checked visually 
every day to collect larvae within 24 h post-hatching. When 
present, larvae were collected by pouring aquarium seawater on a 
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Figure 1. Morphological variables. All measurements were taken on larvae lying on their right side. A: maximum shell length (L) (in urn); B: 
maximum shell height (h) (in um); C: projected surface area of the left side (in mm 2 ); D: deformity index, a is the angle of the abnormality (in degrees) 
and d, the depth of the abnormality (in |im). 
doi:10.1371/journal.pone.0093021.g001 



200 (xm mesh sieve, rinsed with seawater and preserved in 96% 
ethanol. 

Pools of larvae from adults acclimated to the different pC0 2 
levels since January were collected from the different pC0 2 
conditions at the temperature level of 19°C between 8 and 24 June 
2011. Only samples with enough intact larvae were used. Thus 
two viable samples per pCQ 2 condition were studied. 



Shell mineralization 

For each pC0 2 treatment, 5 to 8 larvae were randomly chosen 
among the previous subsamples and observed under polarized 
light to determine birefringence patterns with an Olympus 
dissecting microscope equipped with polarizing filters. All polar- 
ized images were acquired with an Olympus camera at xlOO 
magnification with 40 ms light exposition. Birefringence under 



Morphological variables 

Morphological measurements were performed on a random 
subsample of 40 larvae when possible or at least 20 larvae from 
each of the 6 larval pools. Larvae with unbroken shells were 
isolated in sterile, flat-bottom, 96-well plates and preserved in pure 
glycerol as described in Auzoux-Bordenave et al. [68] . Each larva 
was placed on its right side and photographed under light 
microscopy using an Olympus Camedia C-7070 camera attached 
to an Olympus SZX 12 dissecting microscope. Pictures were taken 
without autofocus at x90 magnification. Maximum length, height 
and projected surface area of the left side (Figures 1A, B, and C, 
respectively) were measured by analyzing images with ImageJ 
software [69], after calibration with a stage micrometer. 

Abnormalities 

In each subsample used for morphological measurements, 
veliger larvae with abnormal shells were counted and the 
percentage of abnormal larvae was estimated per pG0 2 treatment. 
To be considered as a shell abnormality and not as a broken shell, 
deformities had to be devoid of fracture lines. A "deformity index" 
was calculated to quantify the intensity of the shell deformity. It 
was defined as the ratio between the angle formed by the 
abnormality and its "depth", which is the distance between the 
theoretical curve of the shell and the forest point (extreme point) of 
the deformity (Figure ID). 



1000 




390 



750 
pC02 (Matm) 



1400 



~J Length (urn) 
3 Height (urn) 
■ Area (urn 2 10" 3 ) 



Figure 2. Morphological variables. Mean shell length, height and 
left surface area (± SE) in the different pC0 2 treatments. Different 
letters above bars indicate significant differences between treatments 
(p<0.05, Dunn post-hoc test), n = 51 to 92. 
doi:10.1 371/journal.pone.0093021 .g002 
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Table 2. Effect of pC0 2 on morphological variables, abnormality indices and intensity of birefringence. 





Kruskal-Wallis Test 




df 


H 


P 


Length 


2 


37.353 


< 0.001 


Height 


2 


16.235 


< 0.001 


Surface area 


2 


30.106 


< 0.001 


Abnormality index 


2 


6.046 


0.049 


Intensity of birefringence 


2 


14.562 


< 0.001 



Legend: Summary of the non-parametric Kruskal-Wallis tests testing the effect of p€0 2 on each morphological variable, abnormality index and birefringence intensity. 
doi:1 0.1 371/joumal.pone.0093021 .t002 



polarized light is due to the mineral phase composing the shell 
[40,68,70,71]. In the absence of mineralized structures, there is no 
birefringence and the picture looks totally black. Under identical 
light conditions, areas appearing more birefringent contain a 
much larger proportion of crystalline calcium carbonate [70,72]. 
The intensity of birefringence of each shell was used as a proxy for 
mineralization level for the three pC0 2 treatments. It was 
quantified from pictures by using ImageJ software [69]. Pictures 
of polarized shells were first transformed into grayscale images. A 
mean gray value (in pixels) was determined for each birefringent 
zone. All birefringent zones of the shell were compiled to obtain a 
global mean gray value, giving the intensity of the birefringence of 
the whole shell. 

Statistics 

All statistical analyses were performed using the free software R 
2.15.0 version [73]. Normality and homoscedasticity of the data 
were first checked using Shapiro and Levene tests, respectively. 
Due to the non-normality and heterogeneity of variance, the 
influence of pC0 2 on morphological variables, deformity indices 
and birefringence intensity was analyzed using the non-parametric 
Kruskal-Wallis test followed by the Dunn post-hoc test [74]. A 
Chi-squared (x 2 ) test followed by G-tests (likelihood-ratio test) [75] 



were used to compare percentages of anomaly between the three 
pC0 2 conditions. 

Results 

Morphological variables 

pC0 2 significantly affected length, height and surface area of 
the hatched larvae (Figure 2, Table 2). These morphological 
variables are related to each other and were generally influenced 
in the same way by pC0 2 . Length and height were the highest at 
390 itatm and significantly decreased with increased pC0 2 . 
Larvae collected at 750 and 1400 u,atm^C0 2 showed a decrease 
of 5.3% and 5.9% in length, respectively, and 2.6% and 4.5% in 
height, respectively, compared to control larvae (390 ixatm). 
Similarly, the greatest shell surface area was observed at 390 u.atm 
pC0 2 , but then significandy decreased by 6.2% and 1 1.2% at 750 
and 1400 itatm pG0 2 , respectively. 

Abnormalities 

Abnormalities in larvae were observed as notches located close 
to the shell aperture (Figure 3). The percentage of abnormal larvae 
increased with increased pG0 2 and ranged from 6.7 to 26.5% 
(Figure 4; % 2 test, p<0.05). Abnormalities were 1.5- and 4-fold at 
pC0 2 levels of 750 and 1400 ixatm, respectively, than at 




Figure 3. Shell abnormalities. Different intensities of shell abnormalities observed among samples. A, B, C, and D show whole larvae whereas E, F, 

G, and H show the detail of their respective abnormalities. 

doi:10.1371/journal.pone.0093021.g003 
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Figure 4. Percentage of abnormalities in the different pC0 2 
conditions. Different letters above bars indicate significant differences 
between treatments (G-test, p<0.05). 
doi:1 0.1 371 /journal.pone.0093021 .g004 

390 u.atm. Furthermore, different intensities of shell abnormalities 
were observed with variation in notch acuteness. The deformity 
index varied between 0.03 (390 |iatm pC0 2 ) and 0.17 (1400 natm 
pC0 2 ). Although the Kruskal-Wallis test showed a marginally 
significant pC0 2 effect (p = 0.049; Table 2), the pairwise Dunn 
post-hoc test did not detect significant differences between the 
three pC0 2 treatments (p > 0.05). 

Among the abnormal larvae observed under polarized micros- 
copy (see below), some showed abnormalities which appeared less 
birefringent, and even not mineralized, as revealed by the lack of 
birefringence in these parts of the shell (Figure 5). 

Shell mineralization 

Pictures taken under polarized light (Figure 6) suggest that the 
intensity of the birefringence decreased with increased in pC0 2 . 
The measure of birefringence intensity using the mean gray values 
estimated for each shell clearly confirmed this relationship (Table 



2). Mineralization was greatest at 390 u.atm pC0 2 , intermediate at 
750 uatm pC0 2 and lowest at 1400 u.atm pC0 2 (Figure 6). 

Discussion 

In our study, the effects of elevated pC0 2 were integrated from 
embryo formation and throughout embryonic stages up until the 
release of veliger larvae. The integration of/>C0 2 effects across all 
developmental stages, from fertilization to settlement and beyond, 
is particularly instructive because early life stages may respond to 
environmental stressors in a different way than adults. Acute 
exposures of early life stages have shown various effects in growth 
or feeding performances [76,77], but these may not represent field 
conditions. Results obtained from acclimation to high pC0 2 across 
different life stages often differ from those arising from acute 
exposure of a given stage. For example, it has been shown that 
oyster D-veliger larvae grown from eggs fertilized at elevated 
pC0 2 were more drastically affected than those first produced at 
ambient pG0 2 and then reared later (embryo stages) at high/;C0 2 
[78]. Keeping the parents under different pC0 2 conditions before 
fertilization and until larval release allowed embryos develop 
entirely under a given level of stress. To our knowledge, only 
Dupont et al. [79] on sea urchins, Parker et al. [80] on mollusks 
and Vehmaa et al. [81] on copepods acclimated adults to high 
pC0 2 during reproductive conditioning before studying larvae in 
the same pC0 2 conditions. 

pC0 2 effects were first investigated on the shell morphology of 
the hatched larvae. The size (length, height and surface area) of 
the released larvae decreased with increased pC0 2 . Length and 
height were not significantly different between the pC0 2 levels of 
750 and 1400 u.atm, whereas the shell surface area progressively 
decreased with increased pC0 2 from 390 to 1400 u.atm. Although 
in some rare cases, elevated pC0 2 does not affect larval 
morphology and growth [82,83], a correlation between high 
pC0 2 and smaller size has been demonstrated in most bivalve and 
gastropod species studied to date (see review in Table 3), with 
pC0 2 sometimes affecting the shape of the larval shell [41]. As 
observed here for length and height (ca. -5% in the two high 
pC0 2 treatments), shell lengths of Crassostrea gigas veliger larvae are 
smaller under elevated pC0 2 , but are similar in conditions of pH 
lowered by 0.4 and 0.7 pH units, a range corresponding to our pH 
values [32]. Smaller size and delayed shell growth can be 




4 X 



B 



Figure 5. Shell abnormalities under normal (A, C, and E) and polarized (B, D, and F) light. A-B are pictures of a normal larva presenting the 
characteristic black cross of birefringence. C-D and E-F show examples of abnormalities observed among the samples, with arrows indicating less 
calcified zones. 

doi:1 0.1 371 /journal.pone.0093021 .g005 
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Figure 6. Birefringence intensity. Mean gray value (in pixel) of the shells under polarized light in the different pC0 2 treatments (± SE). Bars with 
different letters are significantly different (Dunn post-hoc test, p<0.05,), n = 5 to 8 larvae per condition. Pictures taken under polarized light (below 
the graph) show birefringent patterns in the different pC0 2 treatments. 
doi:10.1371/journal.pone.0093021.g006 



attributed to problems in shell deposition, delayed periostracum 
formation and/or increased shell dissolution, as hypothesized by 
Watson et al. [84] . 

Such processes may lead to developmental abnormalities and to 
an increase in their frequency under elevated />C0 2 . Some (7%) C. 
fomicata larvae in the control group (390 uatm /)C0 2 ) showed mild 
shell abnormalities in the form of a notch close to the aperture. 
The frequency of this abnormality increased under high pC0 2 , 
being 1.5-fold more frequent at 750 uatm pG0 2 and reaching 
26% at 1400 uatm /;C0 2 . The intensity of shell abnormality, 
estimated using the deformity index, did not vary significantly with 
increasing pC0 2 , although more pronounced shell deformities 
were detected at the highest pC0 2 condition (1400 uatm). The 
occurrence of abnormal shells is a common response in mollusk 
larvae exposed to elevated pC0 2 . In bivalves for example, 
abnormalities can occur as shell hinge and edge deformities 
[85] , irregular-shaped shells [40] or protruding mantles [86] . The 
frequency of abnormalities can reach 40% of shell deformities in 
Pecten maximus larvae reared at 1250 uatm pC0 2 [85], and even 
90% of abnormal D-veligers in Saccostrea glomerata at 1000 uatm 
pC0 2 [31]. In gastropods, larval shells are considered abnormal 
when shells are too small to fully cover the soft body [34] or when 
dissolution zones are observed at the edge of the aragonitic larval 
shell [33], with frequencies of abnormality ranging from 20% in 
Haliotis discus hannai at 1650 u.atm pC0 2 [33] to 40% in Haliotis 
kamtschatkana at 800 uatm pC0 2 [36]. At extremely high/)C0 2 (> 
1 700 uatm), some abalone larvae are even unable to precipitate a 
calcareous shell [36,44]. 



Such abnormalities may be due to different processes: (i) the 
production of amorphous CaC0 3 may be affected by damage to 
embryonic ectodermic cells and/or (ii) seawater corrosion may 
induce shell dissolution, affecting the strength and calcification of 
some parts of the shell [32]. Here, the mineralization level of larval 
shells was investigated at each pC0 2 level by observing the veliger 
aragonitic shell under polarized light [70] . The characteristic dark 
cross observed in each larval shell indicated a radial arrangement 
of aragonite crystals [72] and did not have been considered as 
non-crystalline zones. The intensity of birefringence was used as a 
proxy for mineralization because increases in birefringence reflect 
increases in crystalline structure and calcification of the shell. 
Observed under polarized light, abnormalities appeared less 
birefringent than the rest of the shell, suggesting that deformities 
were likely less calcified as proposed by Barros et al. [32] 

The birefringence intensity of the larval shells decreased with 
increased pC0 2 , and was significantly lower at 1400 u.atm^C0 2 . 
This drop in birefringence revealed a decrease in calcification, 
which may be related to a less mineralized matrix [87], or more 
likely to a reduction in shell thickness [72]. Our data did not allow 
us to discriminate between these two possibilities, but previous 
studies have already reported a decrease in shell thickness under 
high pC0 2 in bivalve larvae. For example, using scanning electron 
microscopy measurements, Gazeau et al. [39] showed a decrease 
in thickness of 12% in Mytilus edulis larvae at 745 u.atm pC0 2 . 
Talmage and Gobler [43] report a decrease in thickness of 
Mercenaria mercenaria (-43%) and Argopecten irradians (- 47.5%) larval 
shells after 17 days of development at 1500 uatm pC0 2 , which 
was associated with an impact on the integrity and the 
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Table 3. Literature review of reports of high pC0 2 effects on morphological variables in mollusk larvae. 



pC0 2 conditions % decrease due to 



Reproductive mode 


Species (veliger stage) 


Measured parameter 


pH conditions 


(liatm) 


pco 2 


Study 


Broadcast spawner 


Crassostrea gigas 


shel 


length 


PHnbs 7.4 


2268 


decrease 


[40] 


Broadcast spawner 


Crassostrea gigas 


shel 


length 


pH NBS 7.8 


1000 


16% 


[78] 


Broadcast spawner 


Crassostrea gigas 


shel 


length and height pH NBS 7.76 - 7.37 


1386 - 3573 


10.6% 


[32] 


Broadcast spawner 


Crassostrea gigas 


shel 


area 


pH NBS 7.7 - 7.4 


1497 - 2386 


18.7 - 29% 


[101] 


Broadcast spawner 


Crassostrea virginica 


shel 


length 


pH T 7.84 - 7.49 


<= 650 - 1500 


16.7% 


[38] 


Broadcast spawner 


Crassostrea virginica 


shel 


area 


pH NBS 7.76 


840 


16% 


[94] 


Broadcast spawner 


Saccostrea glomerata 


shel 


length 




1000 


22% 


[31] 


Broadcast spawner 


Saccostrea glomerata 


shel 


length 


pH NBS 7.8 


1000 


34% 


[78] 


Broadcast spawner 


Saccostrea glomerata 


shel 


length 


pH NBS 7.9 


856 


31.6 - 1.3% 


[80] 


Broadcast spawner 


Saccostrea glomerata 


shel 


length 


pH NBS 7.8 - 7.6 


508.8 - 775.6 


8.7 - 6.3% 


[84] 


Broadcast spawner 


Saccostrea glomerata 


shel 


height 


pH NBS 7.8 - 7.6 


508.8 - 775.6 


7.5 - 5.1% 


[84] 


Broadcast spawner 


Mytilus californianus 


shel 


area 


pH NBS 7.75 


970 


5 - 7% 


[37] 


Broadcast spawner 


Mytilus edulis 


shel 


length 


pH NBS 7.8 


1200 


4.5 - 6% 


[39] 


Broadcast spawner 


Mytilus edulis 


shel 


thickness 


pH NBS 7.8 


1200 


12% 


[39] 


Broadcast spawner 


Mytilus edulis 


shel 


area 


pH NBS 7.6 


1388 - 1493 


7-8% 


[16] 


Broadcast spawner 


Mytilus galloprovincialis 


shel 


length and height pH NBS 7.4 


2000 


26 - 20% 


[86] 


Broadcast spawner 


Argopecten irradians 


shel 


length 


pH T 7.84 - 7.49 


= 650 - 1500 


50% 


[38] 


Broadcast spawner 


Argopecten irradians 


shel 


diameter 


pH T (?) 7.5 


1500 


30.7% 


[43] 


Broadcast spawner 


Argopecten irradians 


shel 


thickness 


pH T (?) 7.5 


1500 


43% 


[43] 


Broadcast spawner 


Argopecten irradians 


shel 


length 


pHsw 7.39 


1987 


18.9 - 7.5% 


[102] 


Broadcast spawner 


Pecten maximus 


shel 


length 


pH NBS (?) 7.51 


1627 


10% 


[85] 


Broadcast spawner 


Pecten maximus 


shel 


height 


pH NBS (?) 7.51 


1627 


5% 


[85] 


Broadcast spawner 


Macoma balthica 


shel 


length 


pH NBS 7.7 - 7.2 


1 700 - 4400 


4.3 - 8.5% 


[103] 


Broadcast spawner 


Macoma balthica 


shel 


length 


pH NBS 7.8 - 7.5 




16.2 - 16.9% 


[42] 


Broadcast spawner 


Mercenaria mercenaria 


shel 


length 


pH T 7.84 - 7.49 


»» 650 - 1500 


22.7% 


[38] 


Broadcast spawner 


Mercenaria mercenaria 


shel 


diameter 


pH T (?) 7.5 


1500 


25.5% 


[43] 


Broadcast spawner 


Mercenaria mercenaria 


shel 


thickness 


pH T (?) 7.5 


1500 


43.0% 


[43] 


Broadcast spawner 


Haliotis discus hannai 


shel 


length 


pH T (?) 7.71 


1050 


2.50% 


[33] 


Broadcast spawner 


Haliotis discus hannai 


shel 


length 


pH NBS 7.6 - 7.3 


1 224.6 - 2543.6 


decrease 


[34] 


Broadcast spawner 


Haliotis kamtschatkana 


shel 


length 


pH NBS (?) 8.07 


800 


5% 


[36] 


Brooding in pallial cavity 


Ostrea chilensis 


shel 


thickness 


pH NBS (?) 7.7 - 7.0 




no change 


[50] 


Brooding in pallial cavity 


Ostrea lurida 


shel 


growth 


pH T 7.76 


1000 


5 - 14% 


[104] 


Egg masses 


Littorina obtusata 


shel 


length 


pH NBS (?) 7.6 




10% 


[41] 


Egg masses 


Stylocheilus striatus 


shel 


area 


pH NBS (?) 7.6 




24 - 36% 


[52] 


Egg 


Sepia officinalis 


tota 


weight 


pHT 7.84 - 7.60 


750 - 1430 


no change 


[53] 


Encapsulation + brooding Crepipatella dilatata 


shel 


thickness 


pH NBS (?) 6 




30% 


[51] 


Encapsulation + brooding Crepidula fornicata 


shel 


length 


pH T 7.8 - 7.6 


750 - 1400 


5.3 - 10.7% 


Present study 


Encapsulation + brooding Crepidula fornicata 


shel 


height 


pH T 7.8 - 7.6 


750 - 1400 


2.6 - 13.1% 


Present study 


Encapsulation + brooding Crepidula fornicata 


shel 


area 


pH T 7.8 - 7.6 


750 - 1400 


6.2 - 19.8% 


Present study 



PHn B s : pH on the NBS scale; pH T : pH on the total scale; 
doi:1 0.1 371/journal.pone.0093021 .t003 



pH sw : pH on the seawater 



scale. 



connectedness of the hinge structure. A decrease of 5.7% in shell 
thickness of brooded larvae of the oyster Ostrea chilensis has also 
been observed following a decrease in pH (down to 6.56) within 
the mother's pallial cavity due to valve closure under salinity stress 
[50]. 

The decrease in larval size and mineralization level of the shell 
may be due to reduced CaC0 3 saturation or hypercapnic 
suppression of metabolic pathways involved in the calcification 
process [9]. Very little is known about the conditions occurring 



during intracapsular development and how acidified seawater can 
affect encapsulated embryos. Previous studies have shown that egg 
capsules of some gastropods, including C. fornicata, are permeable 
to water and ions (e.g. [48,71,88]) and it can be assumed that the 
capsule wall in C. fornicata is almost impermeable to gas because of 
its low 0 2 conductance [89]. Under "normal" conditions, 
embryos of C. fornicata will progressively be exposed to hypoxia 
[89] and hypercapnia via their respiration. This may lead to low 
intracapsular pH, as reported in other gastropod species [90], 
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without altering development [57]. Under elevated pC0 2 , 
diffusion of more protons (FT 1 ") from external seawater to the 
intracapsular medium may alter intracapsular carbonate chemis- 
try, thus enhancing metabolic acidosis. Similar acidosis can be 
observed under low salinity stress. For example, a decrease in pH 
to 6.4 recorded within the pallial cavity of the calyptraeid 
Crepipatetta dilatata [47,51] led to the partial shell decalcification 
of the brooded encapsulated embryos [51]. This decalcification 
may cause the release of some carbonate ions (C0 3 2 ~), which 
could bind to free H + to form bicarbonate ions (HC0 3 ), thus 
buffering the intracapsular acidosis and limiting the drastic pH 
effects on larval metabolism. This potential buffering role, in 
combination with a decrease in the CaC0 3 saturation state, is 
likely to affect shell mineralization and calcification. Such 
processes have been suggested to buffer acidosis resulting from 
anaerobiosis in C. fornicata [91], and may also explain our 
observations. Alternative mechanisms that can decrease the 
intracapsular acidosis, such as the active excretion of H + out of 
the capsule through a proton pump, as shown in the cephalopod 
Sepia officinalis [92], need to be investigated. 

Altogether, our results show that, despite the potential 
protective role provided by encapsulation and brooding, elevated 
seawater pC0 2 affected the shells of the released larvae in C. 
fornicata. Embryos of C. fornicata were affected by high/>C0 2 during 
their intracapsular development. However, the overall low 
abnormality rate and low decrease in size suggested they were 
likely less affected than other mollusk early life stages. The natural 
exposure of embryos to low intracapsular pH as demonstrated in 
cephalopod eggs (pH on the seawater scale of the perivitalline fluid 
of ca. 7.35 at 16°C [53]) and gastropod capsules (pH of the 
intracapsular fluid lower than 7 [90]) could confer to C. fornicata 
larvae some resilience to elevated pC0 2 levels. Indeed, it has been 
shown that bivalves naturally exposed to high pC0 2 conditions in 
their habitat (due to high levels of benthic respiration or to 
seawater naturally enriched in C0 2 ) are less affected by ocean 
acidification than other mollusk species [93,94]. Further studies 
are however needed to determine the pH of the intracapsular fluid 
in C. fornicata, and how it will be affected under future scenarios of 
ocean acidification. 
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